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Abstract 

We show that the relativistic dynamics in a Gaussian random electro- 
magnetic field can be approximated by the relativistic diffusion of Schay 
and Dudley. Lorentz invariant dynamics in the proper time leads to the 
diffusion in the proper time. The dynamics in the laboratory time gives 
the diffusive transport equation corresponding to the Jiittner equilibrium 
at the inverse temperature = mc 2 .The diffusion constant is expressed 
by the field strength correlation function (Kubo's formula). 

1 Introduction 

The attempts to define a relativistic version of diffusion started a long time ago 
(see the reviews [1][2]). In [3] [4] it has been shown that relativistic Markovian 
diffusion cannot be defined in the configuration space. It has been discovered 
by Schay [5] and Dudley [6] that Kramers phase space diffusion admits relativis- 
tic generalization. The generalization is unique under the assumption that the 
generator depends only on the phase space variables (x,p) and that the particle 
mass does not change during the random evolution. It seems that the authors 
of ref. [7] unaware of the work of Schay and Dudley have introduced another def- 
inition of the relativistic diffusion which depends on an additional four-vector 
interpreted as a fluid velocity (without this four- vector the model would not be 
covariant). This approach is developed in [1] [8]. In another approach [9]- [10] the 
authors following methods previously applied to the non-relativistic Langevin 
equation obtain the models of refs.[7] [8] as well as the model of Schay and Dud- 
ley. We have discussed [11] drifts added to the relativistic diffusion of Schay and 
Dudley in such a way that an equilibrium is achieved and the detailed balance 
satisfied. This model has been independently discovered and studied in detail 
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in [12]. The existence of an equilibrium which is the Juttner distribution [13] 
has been also shown in the models of refs.[7][8] [9]. These diverse investigations 
did not show that there is a physical model realizing a particular version of 
the relativistic diffusion. The model of Schay and Dudley is appealing by its 
mathematical elegance (it has a generalization to arbitrary pseudoriemannian 
manifolds [14]). Nevertheless, it has been rejected in [1][8] as unphysical. In this 
paper we show that the relativistic diffusion of Schay and Dudley results from a 
dynamical model of a particle motion in a random electromagnetic field in the 
same way as the non-relativistic Brownian motion [15] can be derived from a 
non-relativistic dynamics in the electromagnetic field. Moreover, we show that 
a particular version of the relativistic Ornstein-Uhlenbeck process suggested in 
[11] and [12] (see also [10]) comes from an interaction with the relativistic ran- 
dom electromagnetic field. 

In non-relativistic mechanics the diffusive dynamics can be derived from the 
dynamics in a random force [16] [17] [18] [19] (these papers describe also the his- 
tory of the problem with proper citations). The random force felt by a tracer 
particle can be understood as the force coming from a chaotic motion of other 
particles. The motion in a random electromagnetic field has been studied by 
physicists for a long time because of its relevance to astrophysics, plasma physics 
and high-energy physics [15] [20] [21] [22] [23] [24]. The notion of a random Liou- 
ville operator has been introduced by Kubo [25] . In his paper it has been shown 
that the Markov approximation to the random motion leads to the diffusion 
(for a rigorous proof of the non-relativistic diffusion see [16] [17]). In view of 
the results obtained for non-relativistic dynamics it is suggestive to repeat the 
calculations in the relativistic case in order to see the appearance of the diffusive 
behaviour in a relativistic motion. 

The plan of this paper is as follows. In sec. 2 we define the relativistic dynam- 
ics of a massive particle in an electromagnetic field in the proper time and in the 
laboratory (coordinate) time. We move to a statistical description of classical 
dynamics in terms of the evolution of functions on the phase space generated by 
the Liouville first order differential operator. Such a formalism is a preparation 
for a unifying description of deterministic and random systems. In sec. 3 we de- 
fine a random Gaussian antisymmetric tensor field which can be considered 
as an electromagnetic field whose randomness comes from random sources. To 
the best of author's knowledge the relativistic random electromagnetic field is 
constructed here for the first time. It could be considered as a regular(soft) 
version of quantum electromagnetic field when the short distance singularities 
of QED are ignored. We discuss an evolution of observables generated by (an 
adjoint of ) a random Lorentz invariant Liouville operator (sec. 4) in the sense of 
Kubo [25] . We represent the evolution in a form suitable for averaging over the 
electromagnetic field. We refer to Kubo [25]- [26] for an argument that it is suffi- 
cient to calculate the expectation value of the square of the Liouville operator in 
order to determine the generator of the diffusion which results from the Markov 
approximation of the random dynamics. The calculation of an expectation value 
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of the square of the Liouvile operator generating the proper time dynamics is 
performed in sec. 5. It comes out that till the second order in the electromag- 
netic field the random dynamics coincides with the diffusive dynamics resulting 
from the diffusion generator of Schay and Dudley. This is the relativistic analog 
of the results derived for the non-relativistic motion in refs.[20] [21] [22]. It shows 
that for a particle motion in a random field the relativistic diffusion plays the 
same role as the classical Brownian motion. 

We perform the same calculations in the laboratory time (sec. 6). Now, we 
do not have the explicit Lorentz invariance of the Liouville operator. An average 
of the square of the Liouville operator over the electromagnetic field gives a sur- 
prise. There appears a first order Lorentz non-invariant term. We introduced 
such a drift in [11] as the friction determined by the Jiittner equilibrium dis- 
tribution [13] through the detailed balance condition. The Jiittner distribution 
which comes from the random dynamics corresponds to the inverse temperature 
= mc 2 (m is the particle's mass). In the controversy concerning the rela- 
tivistic analog of the Ornstein-Uhlenbeck process as discussed in [1][2][7] [8] the 
interaction with a random electromagnetic field selects the process determined 
by the detailed balance from the diffusion of Schay and Dudley. 



2 Relativistic dynamics 

The dynamics of a particle with mass m in an external electromagnetic field 
neglecting its own electromagnetic field (produced by the current of a moving 
particle) is described by the explicitly Lorentz covariant four-vector equations 
[27] 

dr^ 1 

=- = — (1) 

dr mc 



mc^ - F^p", (2) 



The vectors are defined in the Minkowski space-time with the metric tensor 
r\^ v = (1, — 1, — 1, — 1), the Greek index fi = 0,1,2,3. We shall also use Latin 
indices j = 1,2,3 . It follows from eq.(2) that 

^" Pll p v )=0. (3) 

Hence, 

p 2 =^' J p li p u = m 2 c 2 , (4) 
From eqs.(l) and (4) it follows that r is the proper time 

dr 2 = dx^dx^. (5) 
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We can eliminate r from eqs.(l)-(2) in favor of x° (we call x° the laboratory 
time). Then, eqs.(l)-(2) read 

dxk 1 k ,n\ 

S> = 7,f' (6) 

k = 1,2,3. 

We move from a study of individual trajectories to functions W (observables) 
on the phase space (x,p). Then, W evolves as 

d rW=—^-F jv ^^j. (8) 
mc ox 1 * mc op 3 

There is no derivative over p° in eq.(8) as p° is expressed by p. Knowing the 
evolution of W we can obtain the trajectory (x(-),p(-)). It is useful to apply a 
statistical formulation of classical dynamics. We begin with a probability density 
fi of the initial points (x,p) in the phase space and determine its evolution 
in time r . We define an expectation value of the observable If in a state 
(probability distribution) f2 as 

n(w) = (ft,w) - J dxdpnw (9) 

Then, the adjoint evolution is 

(n T ,w) = (n,w T ). (io) 

The initial probability distribution ft in the phase space becomes a probability 
distribution of trajectories (#(■), p(-)) satisfying the equation 

d T n = -i£M + F jv £-gi. (ii) 

The evolution of an observable V m the laboratory time is determined by eqs. 
(6)-(7) 

^=P^ 1 Vx^ + ^ P > o 1 0- (12) 

By means of eq.(10) we can again define an evolution of the probability distri- 
bution. Such a statistical description is necessary when we consider the random 
dynamics because in such a case an individual trajectory is not an observable. 
Only averages over trajectories have a physical meaning. 
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3 Random electromagnetic fields 

We consider an antisymmetric tensor field F which is covariant with respect to 
the Poincare group. We assume that an average (•) over F is defined which pre- 
serves the Poincare symmetry. This means that the two-point function defined 

by 

{F pv (x)F ap (x')) = G lxvwp {x - x'). (13) 

is a tensor depending only on x ~ x 1 with the following symmetry properties: 
i)G>„ ;<jp is symmetric under the exchange of indices (pv; x) and (ap; x') and 

antisymmetric under the exchange p — > v and a — > p 

We impose the first set of Maxwell equations (Bianchi identities) on the 

tensor field F 

d^PF ap = 0. (14) 
which imply an equation for the covariance (13) 

In Fourier transforms (F(—k) = F(k) because F is real) eq.(13) reads 

{F pv {k)F ap {k')) = G^, ap {k)5{k - k'), (16) 

where G pu - a p (k) is a tensor which must be constructed from the vector k p and 
the fundamental four-dimensional tensors rj pp and e pupa . Hence, the tensor 
G pv -,ap with the symmetries i) has the form 

Gfiv- i( Tp{k) = di{k){T] p(7 k L ,kp Tjppk v k a + i~iv P k p k a f]vak^k^) (~\'7\ 
+a 2 (k)(ri pa ri lJ p - ri PP ri va ) + a 3 (k)e^ vap , 

where dj are scalars with respect to the Lorentz group. Eq.(15) after the Fourier 
transform reads 

k a e a ^G p ^p(k) = 0. (18) 

It follows from eqs.(17) and (18) that a\ can be an arbitrary function of k 2 and 
a,2 = 0. From eq.(18) we obtain an equation for a% 

k a e a ^e pva pa 3 (k)=0. (19) 

Hence, a,3(k) = if k ^ (i.e. 03 = KS(k)). In the configuration space the 
admissible 03 term would lead to a constant term (K is a constant) 

^ jiv;ap = K- e HV(jp- (20) 

So far we did not make any assumption concerning the question where the 
average value (13) comes from. As an example, it could come from a time 
average. It may result from a certain probability distribution of the currents 

d»F pv = J v (21) 
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or the electromagnetic field becomes random as a result of a passage through a 
random medium. In the latter cases there exists a probability measure deter- 
mining the distribution of F^ v . It is known [28] that a necessary condition for 
an existence of such a measure is 

(( J dxF^{x)r(x)) 2 ) - J dxdx' r{x)ff{x')G llv , ap {x -x')>0 (22) 

for arbitrary test-functions /'"' (satisfying some regularity and growth condi- 
tions [28]). The condition (22) is also sufficient for an existence of a random 
Gaussian field F. 

The positivity condition (22) excludes the e-term (20) because the quadratic 
form 

J dxdx'.r v (x).rP(x')e^ rTp 

cannot be positive for arbitrary test-functions It follows that only the a\ 
term, which is a function of fc 2 , remains in eq.(17). Hence, 

(F^(x)F ap (x')) = -D^, ap G(x - x'), (23) 

where 

= -Vnadvdp + ri w d v d a - i] vp d a dp + r) va dp,d p (24) 

and G is a function depending only on (x — x') 2 . 

There remains to explore the positivity condition (22). We show that under 
the assumptions i)-ii) there exists a random antisymmetric tensor field with the 
two-point correlation function (23) if and only if the Fourier transform G of the 
function G in eq.(23) satisfies the condition 

iii) 

G(k) > (25) 

and G(k)=0 if k 2 < 0. 

We prove that the inequality (22) follows from the condition iii). With the 
two-point correlation function (23) eq.(22) reads 

/ dxdx' f^{x)rP{x')Gp V , r7p {x - x 1 ) 

= - / dk~3^ v (-k)(j) p<J k v k p - r) pp k u k a + r]„pkp,k a - r] VC rkpk p )!' jp G(k) (26) 
= fdkjHkjgi(k)G(k) - / dfcA j / 0j 'fei/°'(*)G(fc), 
where 

9 j =k f 0j + kj l i. (27) 
Expressing kjf ^ by kjg^ we obtain 

/ dxdx'f^{x)r_P_{x')G^ p {x - x') = J dkG(kjg^g3 - J dkG{k)k^ 2 \g^kj\ 2 
>JdkG(k)kv 2 gigi(k 2 -k 2 ) > 0, 

(28) 
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if (25) and the condition iii) are satisfied. 
Remarks: 

f .Note that there is a function M vap such that 

{d»F^(x)F ap {x')) = d^M^ p (x - x'). 

We treat this consequence of eq.(23) as a weak form of the second set of Maxwell 
equations. The condition that the electromagnetic field has no source (J = 0) 
can now be imposed by the requirement that M satisfies the wave equation . 
2. If the source-less Maxwell equations (21) are satisfied then 

G(k) = S(k 2 ) (29) 

(eq.(29), with ko > 0, is satisfied for the vacuum two-point function in quantum 
field theory of the free electromagnetic field [29]). However, the two-point func- 
tion (23) determined by its Fourier transform (29) is singular. The singularity 
would appear in the diffusion equation as a singularity of the diffusion coeffi- 
cients. We do not impose J = in eq.(21) in order to work with more regular 
electromagnetic fields. 

It follows from cqs.(23)-(24) that 

D pv\crpG = T\n<jGvp ^]ppGi/a + T\vpGap TJi/aG^p (30) 

with 

G lil/ (x) = dfj,d u G = rinvhi(x) + x AJ x„/i 2 (x). (31) 

Here, G(x) = g(u) with u = x 2 , hi(x) = 2g'(u) and h 2 (x) = Ag"(u). The 
functions hj(x) depend only on x 2 = x^x^. 

If in a model of an electromagnetic field the scalar G is a function of some 
other vectors (in addition to x) then G is a function of scalars built from these 
vectors. In such a case the formula (31) contains more terms. As a conse- 
quence, the diffusion discussed in subsequent sections would also depend on the 
additional vectors. 



4 Random evolution 

In this section we discuss a solution of the evolution equations (8) and (11). We 
are going to represent the solution in such a form that an averaging over the 
electromagnetic field can be performed. The differential equations (8) and (11) 
for the evolution of functions W on the phase space are of the form 

8 S W S = (X + Y)W S , (32) 

where Y is a random first order differential operator and X is the free evolution. 
For the evolution (8) in the proper time 

X = — jrdp. (33) 
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Let 

Y(s) = c^(-sX)Ycxp(sX). (34) 
Then, the solution of eq.(32) can be expressed as 

W t = exp(tX)W{, (35) 

where 

d s W ! s = Y{s)W I s , (36) 

(this is the interaction picture well-known from quantum mechanics). We can 
solve eq.(36) by iteration. The iteration till the second order reads 

W} = Wo + Si dsY{s)W I s + J* ds j; ds'Y(s)Y(s')W s I , 

= W + J* dsY{s)Wl + § /„' ds J S ds>(Y(s)Y(s>) + Y(s')Y(s) + [Y(s), Y(s')})W^. 

(37) 

We can expand the solution (37) in an infinite series. We assume that the 
initial condition W for eq.(32) is independent of the electromagnetic field. Then, 
the expectation values of the products of Y(s) appearing in the series can be 
calculated because the correlation functions of the Gaussian random held F 
are determined by the two-point function (23). In the form of a path-ordered 
integral the solution has the form 

W/ = T(exp (J dsY(s)^W . (38) 

The diffusion limit applied in sees. 5-6 relies on the approximation [25] [26] 

(T(eMf dsY(s)))) =cxp(i J* ds f ds'((Y(s)Y(s') + Y(s')Y(s)))). 

(39) 

If [Y(s),Y(s')] = and Y is a linear function of Gaussian variables then 
eq.(39) is exact. In our case the operator Y(s) (generating the Liouville dy- 
namics) is an antihermition operator in L 2 (^dx). Then, the anticommutator 
Y(s)Y(s') +Y(s')Y(s) is a hcrmition operator whereas [Y(s),Y(s')] is antiher- 
mition. We neglect the commutator in eqs.(37) and (39) corresponding to the 
unitary Liouville evolution and restrict ourselves to the pure diffusion. 

5 The expectation value of the proper time evo- 
lution 

The observable time evolution results from an averaging over fluctuations of the 
electromagnetic field. The fluctuations of the electromagnetic field make the 
dynamics chaotic. The particles' trajectories become random. In this section 
we apply the approach of Kubo [25] [26] in order to obtain a diffusion as a Markov 
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approximation to random dynamics. The result involves an approximation of 
the Liouville evolution of functions on the phase space by an evolution generated 
by second order differential operators. In this formulation it is sufficient to 
calculate an expectation value of the square of the Liouville operator in order to 
determine the generator of the diffusion. We consider the proper time evolution 
(l)-(2) first. In general, we could split the Liouville operator Y tot = Y ex + Y, 
where Y ex is the Liouville operator corresponding to an external deterministic 
electromagnetic field and Y describes the random part. The Y ex part could be 
added to the final result. We restrict our discussion to the random Liouville 
operator (34). Taking as Y the second term in eq.(8) and X from eq.(33) we 
calculate Y(s) (34) 

Y(s) = F jv (x-—p)p v ^ 1 . (40) 
mc op> 

We apply the covariance (30)-(31) (with x M ~ ^p^ from eq.(40)) to calculate 
the expectation value appearing in eq.(39) 

I So ds So ds'((Y(s)Y(s>) + Y(s')Y(s))) = (mc)-* f Q ds f° ds' 
(vjl(Vv P H i(s - s') + m- 2 cr 2 p v p p H(s - s')) 

-r) jp (r) u iHi(s - s') + m- 2 c - 2 p l ,piH(s - s')) (41) 
+V<y P (VjiHi(s - s') + m- 2 c- 2 p : jp l H(s - s')) 

-M^H^s - s>) + m- 2 c- 2 p jPp H( S - s'^p^pP^, 

where 

H 1 (s-s') = h 1 {{s-s') 2 ) = 2g'(u) 

and 

H(s - s') = (s- s') 2 h 2 ((s - s') 2 ) = Aug"(u). 

where u = (s — s') 2 , hj(x) (defined after eq.(31)in terms of G(x 2 ) = g(Vx 2 )) 
being functions of x 2 depend only on s — s' because after a calculation of the 
two-point function (13) of F(x — -^p) and F(x — -^.p) we obtain 

h j (((s-s')p(mc)-'f) = h j ((s-s>n 
The explicit calculations in eq.(41) give 

I f r ds J S ds'((Y(s)Y(s') + Y(s')Y(s))) = (mc)- 2 J Q T ds J Q S ds' ^(s - s') + H(s - *')) ^ 

(42) 

where 

AS = (P l m 2 c 2 +p>pl)^^ + 3p 1 ^ (43) 

is the generator of the relativistic diffusion of Schay [5] and Dudley [6] . In the 
approximation (39) (neglecting the non-commutativity) the result of an average 
over a random electromagnetic field is the relativistic diffusion with the time 
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transformed by the function standing in front of in eq.(42). Note that 

is self-adjoint in L 2 (dppQ 1 ) similarly as its non-relativistic analog (generator of 

Krammers diffusion) is self-adjoint in L 2 (dp). 



6 Random evolution in laboratory time 

Eqs.(l)-(2) are equivalent to eqs.(6)-(7). Both sides of eqs.(l)-(2) are four- 
vectors transforming linearly under the Lorentz group. This property does not 
hold true for eqs.(6)-(7). The equivalence between the description in terms 
of the proper time and coordinate time may be lost if we perform non-linear 
transformations such as the averaging over random fields. For this reason we 
calculate the expectation values of the square of the laboratory time Liouville 
generator (12) again. In the laboratory time the free evolution is determined by 

X=p - 1 pV x . (44) 

Then from eqs.(12) and (34) 

Y(s) = F iAJ (x - sp-ip.s^p" 1 ^. (45) 

For the calculation of the expectation values of Y (45) we must make the re- 
placements in formulas of sec. 5 

x -> Po X ps 

and xq — x' = s — s'. So that 

{x-x'f = p^ 2 m 2 c 2 (s-s') 2 (46) 

as the argument of the functions hj in eq.(41). From eq.(46) we can see that 
the transition from the proper time evolution to the laboratory time evolution 
involves p — > PPq 1 ■ This is equivalent to s — > smcp^ 1 in the functions hj when 
we calculate the expectation value (41). Changing the time integration variables 
s — > smcpQ in eq.(41) we obtain (we write xq = ct) 

\ C ds fSds>{(Y(s)Y(s>) + Y(s')Y(s))) 

= j^Po 1 ds f* ds'^rj^H^s - s') + m- 2 c- 2 p„p p H{s - s')) 

-rjj P (VvlHi((s - s') + mr 2 c- 2 p„ Pl H( S - «')) ( 47 ) 
+7]u P (VjlHi((s - s') + m- 2 c- 2 p 3 piH(s - s')) 

-Vvl(Vj P Hi(s - s') + m- 2 c- 2 p 3 p p H(s - s')j)pop" ^"Pq 1 ^r- 
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We have an additional term D\ (coming from p> differentiation of p 1 on the 
rhs of eq.(47)) in comparison to the rhs of eq.(41) 



D 1 = - Po -VV — . (48) 



This term after a contraction with the tensors in eq. (47) gives the result 



§ /* ds /" ds'((Y(s)Y(s') + Y(s>)Y(s))) = ds J Q S ds'^H^s s') + H{ 

(49) 



where the last term in eq.(50) is the friction introduced in [11]. gener- 
ates a diffusion which equilibrates to the Jiittner distribution with the inverse 
temperature P^ 1 — rac 2 . We obtain the Jiittner equilibrium distribution [13] 
= cxp(— (3cpo) from the requirement 



where the adjoint is in L 2 (dp) (see eq.(9)). 

It is surprising that the proper time dynamics (l)-(2) and the laboratory 
time dynamics (6)- (7) after an average over the electromagnetic field give in- 
equivalent results. The probability distribution f2 of the random dynamics in 
the laboratory time tends to the Jiittner equilibrium distribution $ e when time 
goes to infinity. The covariant dynamics of sec. 5 has no limit when the proper 
time or the laboratory time go to infinity. It is easy to comprehend the differ- 
ence from the technical point of view as the Liouville generator in laboratory 
time dynamics has the Pq 1 factor which leads to the extra term (48). However, 
some mysteries remain concerning the question why equivalent dynamics have 
different diffusive limits. 

We assume that the functions Hi and H decay fast for s^s'. We approxi- 
mate the s'-integral appearing in eqs.(42) and (49) as follows 



/ S ds'(2H 1 (s - s') + H(s - s')) ~ f °° ds'(2H 1 (s') + H(s')) = (51) 



where k 2 is the diffusion constant of rcf.fll]. We apply the formulae below 
eq.(41), expressing Hi and H in terms of G{x) to calculate 



Jo 

,where G{x) = g(u) with u = x 2 . 

The Markov approximation must be performed also for higher order terms in 
the expansion (37) in order to justify the formula (39). The approximation has 



with 




(50) 



= 0, 




(52) 
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been discussed first by Kubo [25]-[26] (the expression of the diffusion constant 
(51) and (52) by the correlation function of the forces defining the functions 
hj is known as the Kubo formula). Summarizing the results (42) and (49) we 
have in the Markov approximation 

d r Wl = ^A%W! (53) 
When we define W as in eq.(35) then W satisfies the equation 

d T W T = P^W T + Y A H W r- (54) 

In the laboratory time the corresponding equations read 
and 

^m* w = Y Af>w - (56) 

7 Summary and outlook 

We have been studying a model of particle relativistic dynamics in a random 
electromagnetic field showing diffusive behaviour. The aim was to find which of 
the models of the relativistic diffusion encountered in the literature is realized in 
this physical model. For this reason we have calculated an expectation value of 
the square of the Liouville operator which defines the generator of the diffusion. 
In the case of the dynamics in proper time, as the result of the calculations, 
we obtain the generator of the diffusion of Schay and Dudley; a mathematical 
model discovered a long time ago but disregarded by most physicists. We repeat 
the calculations for the dynamics in the laboratory time. As a result we obtain 
the generator of the proper time diffusion plus a drift term which moves the 
probability distribution to the equilibrium (the Juttner equilibrium distribution 
at /3 _1 = mc 2 ). Comparing with the non-relativistic dynamics in the electro- 
magnetic field we can see that the diffusive dynamics in the proper time can 
be treated as a relativistic analog of the Brownian motion whereas the dynam- 
ics in the laboratory time is analogous to the one of the Ornstein-Uhlenbeck 
process. It follows that the diffusion of Schay and Dudley should be treated 
as a starting point for an investigation of the diffusive motion of relativistic 
systems in the same way as the Brownian motion (in Krammers version) serves 
for a construction of other non-relativistic diffusions in physics [15] as well as in 
mathematics. 

Some extensions of our formulation can be suggested. Applying the Martin- 
Siggia-Rose formalism [30] [18] [19] one could obtain a functional measure (de- 
fined on trajectories) corresponding to the motion in a random electromagnetic 
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field. The exact formula for the functional measure can be applied for fur- 
ther investigations of the random motion beyond the Markov approximation. It 
could also be used for an inquiry of the meaning of the proper time dynamics in 
random relativistic systems and the relation to laboratory time dynamics. The 
equilibration to the Juttner distribution (expected on physical grounds) result- 
ing from the generator (50) may indicate a deeper mathematical mechanism (in 
the path integral formalism) leading to the relativistic version of the Ornstein- 
Uhlenbeck process. In another direction, we can assume a dependence of the 
two-point function (13) on some other tensors (e.g. a vector describing a mov- 
ing frame for theories at finite temperature) . In such a case the formalism does 
not change. We can calculate the dependence of the resulting diffusion (which 
can be treated as a perturbation of the one studied here) on the additional 
variables. Then, astrophysical applications of the relativistic diffusion require a 
general pseudoriemannian metric (as in [14]). The formulation of sees. 2-6 can 
be extended to general relativity. A local Markov approximation leading to the 
diffusion on a manifold is still possible. However, a global approximation of the 
relativistic dynamics on a manifold by the diffusion of ref.[14] may encounter 
some difficulties. 

References 

[1] C. Chevalier and F. Debbasch, 
AIP Conf.Proc.913,42(2007) 

[2] J. Dunkcl and P. Hanggi,Phys.Rep.471, 1(2009) 

[3] J. Lopuszanski, Acta Phys.Pol. 12,87(1953) 

[4] R. Hakim, J.Math.Phys.9, 1805(1968) 

[5] G.SchayPhD thesis, Princeton University,1961 

[6] R.Dudley, Arkiv for Matcmatik,6,241(1965) 

[7] F. Debbasch, K. Mallik and J.P Rivet, Journ.Stat.Phys.88,945(1997) 

[8] C. Chevalier and F. Debbasch, J.Math.Phys.49,043303(2008) 

[9] J.Dunkcl, P. Talkner and P. Hanggi, Phys.Rev.D75,043001(2007) 

[10] J. Dunkcl and P. Hanggi, Phys.Rev.E72,036106(2005) 

[11] Z. Haba, Phys.Rcv.E79,021128(2009) 

[12] J. A. Alcantara and S. Calogcro, Kinetic and Related Models, 4,401(2011) 

[13] F. Juttner, Ann.Phys.(Leipzig)34,856(1911) 



13 



[14] J. Franchi and Y. Le Jan, Commun.Pure and Appl.Math. 60,187(2007) 

[15] S. Chandrasekhar, Rev.Mod.Phys.15, 1(1943) 

[16] H. Kesten and G.C. Papanicolaou, Commun.Math.Phys.78, 19(1980) 

[17] T.Komorowski and L. Ryzhik, Commun.Math.Phys.263,277(2006) 

[18] N. Lebedev, P. Maass and S. Feng, Phys. Rev. Lett, 74,1895(1995) 

[19] L. Golubovic, S. Feng and Fan-An Zeng, 
Phys. Rev. Lett. 67,2115(1991) 

[20] J.R. Jokpii,Astrphys.J.146,480(1966) 

[21] P.A. Sturrock, Phys.Rev.141, 186(1966) 

[22] D.E.Hall and P.A. Sturrock, Phys.Fluids 10,2620(1967) 

[23] W.B. Thompson and J. Hubbard, Rcv.Mod.Phys.32, 714(1960) 

[24] B.Svetitsky, Phys.Rev.D37,2484(1988) 

[25] R. Kubo, J.Math.Phys.4,174(1962) 

[26] R. Kubo, M. Toda and N. Hashitsume, 

Statistical Physics II.Nonequilibrium Statistical Mechanics, Springer, 
Berlin, 1985 

[27] L.D. Landau and E.M. Lifshits, Field Theory, Pergamon Press, 
New York, 1981 

[28] I.M. Gelfand and NY. Vilcnkin, Generalized Functions. Vol. IV, Academic 
Prcss,New York,1964 

[29] S.S. Schweber, An Introduction to Relativistic Quatum Field Theory, 
Row, Peterson and Co, Evanston,1961 

[30] P.C. Martin, E.D. Siggia and H.A. Rose, Phys.Rev.A8,423(1973) 



14 



